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ABSTRACT 

The two leading models for the structure of GRB jets are the uniform jet 
model and the universal structured jet (USJ) model. In the latter, all GRB jets 
are intrinsically identical and the energy per solid angle drops as the inverse 
square of the angle from the jet axis. The simplicity of the USJ model gives it 
a strong predictive power, including a specific prediction for the observed GRB 
distribution as a function of both the redshift z and the viewing angle 9. We show 
that the current sample of GRBs with known z and estimated 9 does not agree 
with the predictions of the USJ model. This can be best seen for a relatively 
narrow range in z, in which the USJ model predicts that most GRBs should be 
near the upper end of the observed range in 9, while in the observed sample most 
GRBs are near the lower end of that range. Since the current sample is very 
inhomogeneous (i.e. involves many different detectors), it should be taken with 
care and cannot be used to rule out the USJ model. Nevertheless, this sample 
strongly disfavors the USJ model. Comparing the prediction for the observed 
GRB distribution both in 9 and in z, with a larger and more homogeneous GRB 
sample, like the one expected from Swift, would either clearly rule out the USJ 
model, or alternatively, provide a strong support for it. The test presented here 
is general, and can be used to test any model that predicts both a luminosity 
function and a luminosity-angle relation. 

Subject headings: gamma-rays: bursts — ISM: jets and outflows — radiation 
mechanisms: nonthermal 
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1. Introduction 

Since the discovery of gamma-ray burst (GRB) afterglows in early 1997, several lines of 
evidence have emerged in support of collimated outflows, or jets. The total energy output in 
7-rays that is inferred from the fluence of bursts with a known redshift z, assuming spherical 
symmetry, E 7jiso , in some cases approaches and in one case (GRB 990123) even exceeds 
M Q c 2 . This is hard to produce in any progenitor model involving a stellar-mass object. A 
jet can significantly reduce the total energy output for a given E lt i so . A more direct line of 
evidence in favor of jets in GRBs (and probably the best evidence so far), is from achromatic 
breaks in the afterglow light curves (Rhoads 1997,1999; Sari, Piran & Halpern 1999). 

However, despite six years of extensive afterglow observations, the structure of the 
relativistic jets that produce GRBs is still unknown. The structure of GRB jets is a very 
fundamental and important property which effects the requirements from the source that 
accelerates and collimates the jets, and has direct bearing on two of the most basic properties 
of any astrophysical radiation source: the rate and the total amount of energy output. Within 
the fireball model (for review see Piran 2000, Meszaros 2002) there are two very different 
jet structures which are compatible with the observations: (i) the uniform (or 'top hat') 
jet model (Rhoads 1997,1999; Panaitescu & Meszaros 1999; Sari, Piran & Halpern 1999; 
Kumar & Panaitescu 2000; Moderski, Sikora & Bulik 2000; Granot et al. 2001,2002), where 
the initial energy per solid angle e and Lorentz factor T are uniform within some finite 
half-opening angle 9j and sharply drop outside of 9j, and (ii) the universal structured jet 
(US J) model (Postnov, Prokhorov & Lipunov 2001; Rossi, Lazzati & Rees 2002; Zhang & 
Meszaros 2002), with a standard jet structure for all GRBs where e oc 9~ 2 (outside of some 
core angle). It is important to note that in the US J model all GRB jets are intrinsically 
identical (both in their angular profile and total energy). Both jet structures can explain 
the observed correlation between i? 7) i S o and the jet break time tj in the optical afterglow 
light curve, tj oc E~- so (Frail et. al. 2001; Bloom et al. 2003). In the US J model this 
determines the jet structure (e oc 9~ 2 ). In the 'top hat' model tj depends mainly on 9j 
(Rhoads 1997,1999; Sari et. al. 1999), while in the US J model it depends mainly (and in a 
similar way) on the viewing angle # b s from the jet axis. 

The simplicity the USJ model gives it a strong predictive power. In a recent paper 
Perna, Sari & Frail (2003; hereafter PSF03) used this feature of the USJ model to predict 
the observed distribution of viewing angles, n(9) = dn/d9 (hereafter we use 9 instead of 
#obs for brevity). They have shown that the current limited sample of 16 bursts with known 
9 and redshift z fits the predicted distribution very well. In this Letter we extend the 
comparison between the USJ model predictions and observations into one more dimension - 
the redshift z. Namely, we use here the two dimensional (2D) distribution n(z, 9) = dn/dzd9 
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and compare it with the known 9 and z of the observed sample. The ID distribution that 
was used by PSF03 for comparison with the data is n{9) = f n(z, 9)dz. However, the known 
redshift of these GRBs is not used in the ID analysis. Comparing the data both for z and 
for 9 with the 2D distribution, n(z,9), reveals that the agreement between the data and 
the ID distribution, n(9), that was found by PSF03 is accidental, and arises because of the 
integration over the z. The 2D data shows a very poor agreement with the model, and the 
hypothesis that the data is drawn from the model is rejected at 99% significance by a 2D K-S 
test. Thus, the agreement between the data and the ID distribution is misleading, and does 
not provide support for (and certainly does not prove) the USJ model. In order to test this 
statistically, while trying to minimize the possible selection effects as a function of z, it is 
most appropriate to compare n(z, 9) with the data over a relatively narrow range in z. Such 
a comparison has another important advantage: it depends only weakly on Rgrb(z) - the 
GRBs rate as a function of z - which is rather poorly known. A drawback of this test is that 
it requires a large number of data points at a given redshift. The narrowest range in z that 
contains a reasonably large number of the current data points (10 points) is 0.8 < z < 1.7. 
In this range the USJ model fails to explain the data with a significance of 99.8%. 

We perform several tests in order to check the robustness of this result, and its sensitivity 
to various selection effects that we can quantify, and find it to be robust and significant. 
Thus, we find that the current data set strongly disfavors the USJ model. Nevertheless, 
it is still premature to draw a definite conclusion, mainly because many different detectors 
were involved in detecting the current sample. This situation is expected to improve in the 
near future with the launch of Swift. Once a homogeneous and large sample of bursts with 
measured z and estimated 9 would be available, applying the 2D test described here would 
result in either a definite rejection or a strong support for the USJ model. We stress that the 
specific test we carry here is relevant only to a universal structured jet, i.e. with a universal 
profile of both e(9) and the initial value of T(^). 1 Finally, although in this Letter we apply 
the 2D test only to the USJ model, it can be easily generalized to any model that predicts 
a luminosity function and a luminosity-angle relation. 



2. Theory 

Below we follow PSF03, and generalize their ID distribution n(9) to the 2D distribution 
n(z, 9). Eq. 5 of PSF03 presents the photon peak luminosity in the energy range 50—300 keV, 



1 if the initial T profile is not universal, then different jets may produce 7-rays within different solid angles. 
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assuming that all GRBs are identical with a differential photon spectral index, a — 1, 

L ph (9,T) = 1.1 x 10 57 f- l 9- 2 photons sec" 1 , (f) 

where T = T sec is an "effective" duration which is given by the ratio of the (isotropic 
equivalent) energy output and peak luminosity (or equivalently (1 + z)~ l times the ratio of 
the fluence and peak flux). In practice, T changes from one burst to another, and we will 
denote its probability distribution by P(T). For a detector with a given limiting flux for 
detection F p h,ii m = -^ph,iim photons cm" 2 sec" 1 and a burst at given 9 and z, we can derive 
the maximal T for which this burst is detected, T max , 

T max = 88 (1 + z)- a D^(z)F ph ] hm (JL^j sec , (2) 

where D 2 &{z) is the co-moving distance in units of 10 28 cm. The total rate of bursts with an 
inferred viewing angle between 9 and 9 + d9 and a redshift between z and z + dz is then 

= * = sin „^S5(£)^ / T "" M p(T)dT , (3) 
v ' ; dzd9 (1 + z) dz J K J ' v ; 

where Rgrb(z) is the GRB rate per unit co moving time per unit comoving volume V(z). Eq. 
3 is similar to Eq. 11 of PSF03 but without integration over z. Thus, Eq. 3 describes the 
2D distribution n(z,9) while Eq. 11 of PSF03 describes the ID distribution n{9) = dn/d9. 

In order to calculate n(z,9) one must assume some Rgrb(z) and P(T). Below we 
consider the Rgrb(z) which PSF03 used as their "standard" model 2 where for z < 10: 

r>ti • (4) 

and ^peak = 2. At z > 10 the rate declines rapidly. PSF03 first considered a delta function 
in T, P(T) = S(T — To), where To is a free parameter which they used in order to get a good 
fit to the observed dn/d9. Using their "standard" Rgr&(z), they found a best fit value of 
T = 8 sec. However, since P{T) can be estimated pretty well from observations, we do not 
take it as a free function. In order to estimate P(T) we used the flux table of the BATSE 4B- 
catalog. We have used the peak fluxes (in photons cm" 2 sec" 1 ), averaged over the 1024 ms 
BATSE trigger, and the fluences from the catalog in the energy range 50 — 300 keV. In 
order to convert the peak fluxes to erg cm" 2 sec" 1 we used a Band spectrum (Band et al. 



2 This model is the Rowan-Robinson (1999) star formation rate with a cutoff at z > 10 as seen in the 
numerical simulation of Gnedin & Ostrikcr (1997) 
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1993) for the energy distribution, with a = -1.0, (3 = -2.0 and E = 100 keV. Finally, T is 
approximated by the ratio of the fruence and the peak flux. 3 According to this estimation 
of T, the total distribution P(T) is consistent with a lognormal distribution, 
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(5) 



with pi = 2.15 and a\ nT = 0.87 (T = 8.6±g 2 sec). Below we first use the delta function for 
P(T) that was used by PSF03, and later we use our Eq. 5. 

Next we consider measurement errors or intrinsic scatter in e9 2 . The measurement 
errors in z are typically negligible. The error in 9, however, may be as large as tens of 
percent, and so is the intrinsic scatter in e9 2 . In order to account for this scatter, and since 
it is more reasonable to assume a Gaussian scatter in In 9 rather than in 9 (both for error and 
intrinsic scatter), we change variables from 9 to ln#, and convolve n(z, In 9) = dn/ dzd\n9 = 
9dn/dzd9 = 9n(z,9) along the In 9 coordinate with a Gaussian of standard deviation <r ln0 , 



n(z, In 9) 



dn _ C iw/2) n(z, ln9') exp [-(ln0 - \n9') 2 /2a 2 n9 ] d(\n9') 



dzd In 9 



J ln(V2) exp [-(ln0 - \n9') 2 /2a 2 ne ] d{\n9>) 



(6) 



Now n(z, In 9) is a rate function which is smoothed along the ln9 dimension by the typical 
scatter due to the bursts intrinsic properties and 9 measurement error. The total scatter 
cannot exceed the measured scatter in e9 2 that was found by Frail et al. (2001). 



3. Results 

First we repeat the ID analysis of PSF03 in 2D for their "standard" Rqrb{z) model 
(Eq. 4), using P(T) = 5(T — T ) with T = 8 sec, for which they obtained the best fit to 
the data. We used the same parameters as PSF03: -F p h,iim = 0.424 photons cm~ 2 sec -1 (The 
threshold for the BATSE trigger on 1024 ms; Mallozzi, Pendleton & Paciesas 1996), a = 1, 
and the same cosmology: VL M = 0.3, Ha = 0.7, and H = 71 km sec -1 Mpc" 1 . 

Fig. la depicts the 2D distribution n(z, ln9), and the circles mark the 16 GRBs of the 
Bloom et al. (2003) sample, which were used by PSF03. This figure is a 2D representation of 



3 This ratio includes the effect of cosmological time dilation, and is therefore a factor of (1 + z) larger than 
the actual value of T, which is measured at the cosmological frame of the GRB. Therefore, we overestimate 
both (T) (by a factor of ~ 2 — 3) and a\ n T (as part of the observed scatter is due to the scatter in (1 + z) 
between different GRBs). Both affects worsen the fit between the US J model and the data. 
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Fig. 1 of PSF03. When integrating our Fig. la over the z dimension we reproduce Fig 1 of 
PSF03 (model 1). Fig. la shows that while the ID distribution ,n(0), provides a good fit to 
the data, the 2D distribution, n(z, 9), does not agree with the data. The 2D K-S test rejects 
the null hypothesis that the data is drawn from the model with a confidence of 99%. The 
reason for the striking difference between the 2D analysis and the ID analysis is that the 
data disagree with the 2D model in two ways that roughly cancel out when integrated over 
redshift: (i) at high z (> 2) there are not enough bursts with low 9 , and (ii) at low z there 
are too many burst with low 9 compared to the number of bursts with high 9 (> 0.2 rad). 
When integrating over the redshift these two shortcomings roughly cancel each other out. 
The fact that the ID distribution of PSF03 peaks at 9 ~ 0.12 rad arises from contribution 
of predicted bursts at high z and low 9, which are not present in the observational sample, 
but are compensated for by the overabundance of bursts at low z and low 9. 

This disagreement with the data, however, cannot be used to draw strong conclusions. 
The reason is that the current sample suffers from numerous selection effects, mainly in 
redshift. The selection effects in z can be minimized by testing the 9 distribution for a given 
z (i.e dn/d9 for a given z). This test has another important advantage: it depends only 
weakly on the poorly known GRB rate -Rgrb(z). The main disadvantage of this test is that 
the size of the data sample is reduced. We therefore take a slice in redshift of 0.8 < z < 1.7, 
which contains 10 of the 16 bursts in the current sample. We use -Rgrb(^) and P(T) from 
Eqs. 4 and 5. We account for a 20% scatter in \n9 by using n(z, In 9) from Eq. 6 with 
<j\ n8 = 0.2. Fig. lb shows the expected distribution n(z,\n9) in this redshift range. Here 
the paucity of bursts with large 9, and overabundance of bursts with small 9, is clear. The 
concentration of bursts at 9 < 0.1 while 9 max (z) ~ 0.25 — 0.4 in this z range contradicts the 
predictions of the US J model - n(9) oc sin#, for 9 < 9 max . The 2D K-S test rejects the model 
with a confidence level of 99.8%. 

In order to check the reliability and robustness of our results, we consider below the 
sensitivity of the results to our assumptions and to the values of the different parameters 
(-Speak, ^ph,iim, etc.). We also consider different observational selection effects and carry 
additional tests to estimate their influence on the results. First, we varied the value of 
Zpeak (Eq. 4), and found that it has almost no influence on the results. We considered also 
the possibility of a larger cr ln . 4 Bloom et al. (2003) obtain a factor of 2.2 scatter in e9 2 , 
implying ai n g < 0.4. Repeating our analysis with a\ n g = 0.4, the US J model is rejected at 
98.7% confidence. 

Next we consider the dependence of our result on the value of -F ph lim . It is important 



4 As discussed above it accounts for both the intrinsic scatter and measurement errors. 
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to note that we consider here a stiff threshold (i.e. constant -F p h,iim) • In reality it is not 
the case (both the detection threshold of the detector and the level of the background 
vary). This effect may result in an underestimate of the number of weak events. One 
way to overcome this obstacle is by choosing a relatively bright threshold which is above 
the detection threshold at any time. Unfortunately, this significantly reduces the size of 
the observed sample, preventing us from applying this method here (it may be applied to 
a future larger sample, e.g. Swift). The only test we can do with the current sample is 
to check the effect of a larger (stiff) threshold. Lower sensitivity (larger F phjlim ) reduces 
Omax(z). In order for the model to accommodate the data concentration at 9 < 0.1 rad and 
0.8 < z < 1.7, Fp h i im needs to be increased by a factor of ~ 5. However, this would imply 
9max(z — 2) m 0.06 rad, which is significantly inconsistent with the observational values of 9 
at this redshift (two bursts with 9 ~ 0.12 and two with 9 ~ 0.22). Another threshold related 
selection effect is the low sensitivity of BATSE in the X-ray. This effect is important because 
of the correlation E p oc e 1 / 2 (Amati et. al 2002), which in our context implies E p oc 9~ l . 
Thus, BATSE is less sensitive to bursts with large 9 (low E p ). Five bursts from our specific 
sample are also used by Amati et al. (2002): four with 9 < 0.1 rad and E p > 400 keV, 
and one (GRB 970508) with 9 = 0.38 rad and an intrinsic E p ~ 150 keV. This sub-sample 
roughly follows the relation E p oc 9~ l , and demonstrates that in the range of 9 where there is 
a deficit of observed bursts, ~ 0.2 — 0.3 rad, the expected intrinsic E p is ~ 200 keV (observed 
E p ~ 100 keV), which is well within the range of BATSE. Therefore, although we cannot 
quantify this effect accurately, it should not strongly affect our sample with limited z range. 

Next we consider the selection effects in 9. A very small 9 implies a very early jet break 
time tj, which may be before the first optical detection and thus result only in an upper limit 
on 9. A large 9, on the other hand, results in a late jet break time tj, which occurs when 
the optical afterglow is too dim for detection (it can be either dimmer than the detection 
threshold or its host galaxy). This would result in only a lower limit on 9. The sample of 
Bloom et al. (2003) shows both effects: in the redshift range 0.8 < z < 1.7 there is one burst 
with an upper limit on 9 and two with lower limits on 9. In order to account for the above 
selection effects, we added the latter three bursts in the most favorable way for the USJ 
model. Namely, we assigned to each of the three bursts the value of 9 within the allowed 
range where n(9, z) assumes its maximal value. Even after adding these three data points 
in this way, the 10+3 data set rejects the USJ model with a confidence level of 99.8%. 

Another selection effect in 9 may result from "dark" bursts. These are bursts with an 
observed X-ray afterglow in which an optical afterglow was not detected despite deep and 
rapid follow-up observations. De Pasquale et al. (2003) argue that a large fraction of the 
"dark" bursts are intrinsically dim, by showing that they have, on average, dimmer X-ray 
afterglows. In the USJ model the intrinsically dim bursts may be interpreted as bursts with 
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very large 9 which are barley detected in 7-rays and X-rays. Since about 50% of the bursts 
with observed X-ray afterglows have no detected optical afterglow, we will make the extreme 
assumption that all of these bursts are intrinsically dim (i.e. have a large 9). We estimate 
the effect of this assumption by adding 13 fictitious points (to the 10 + 3 original data points) 
in the most favorable way for the model. Namely, we generated z according to the model z 
distribution, and then for each point we choose 9 where n(9, z) is maximal (at the given z 
of that point). Even after making these extreme assumptions, the 26 "data" points reject 
the model with confidence level larger than 90%. We conclude that our result, that the US J 
model is incompatible with the current data, is robust and significant. 



4. Discussion 

We have shown that the strong predictive power of the universal structured jet (USJ) 
model enables a determination of the expected distribution of observed GRB rate as a 
function of both redshift z and viewing angle 9. We have compared this predicted 2D 
distribution to current observations, and found a very poor agreement. This is in contrast 
with the result of PSF03 which compared only the observed distribution of 9 to the ID 
prediction of the USJ model (that is obtained from the 2D prediction by integrating over z) 
and found a good agreement. Our analysis shows that this agreement is accidental (resulting 
from the integration over z) and does not support the USJ model in any way. 

However, the poor agreement between the data and the 2D distribution should be taken 
with care, and may not be used to draw definite conclusions. This is since the current GRB 
sample is highly non-homogeneous. It involves many different instruments, and is likely 
effected by various selection effects, which are hard to quantify very accurately. A larger 
and much more homogeneous sample of GRBs with known z and 9 is expected with the 
upcoming launch of Swift, which would enable much stronger and clearer conclusions to be 
drawn from a similar statistical analysis as was done in this work. 5 This would either clearly 
rule out or strongly support the USJ model. Nevertheless, we point out that at least some 
of the selection effects may be overcome by restricting the analysis to a relatively narrow 
range in z. This significantly reduces the redshift selection effects and the uncertainty that 
is introduced by the assumption that has to be made about the poorly known GRB rate 
Rgkb{z)- The main drawback of this second method is that in order to obtain a statistically 
significant sample in a very narrow range in z, many GRBs (many more than in the current 



5 Note that the effect of the variable detection threshold and the E p — e correlation, should be considered 
carefully in the analysis of Swift results as well. 
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sample) are required. Again, such a sample is expected to become available with Swift. 

In the current sample, the narrowest redshift range which still contains enough data 
points for the purpose of statistical analysis is 0.8 < z < 1.7 (it contains 10 out of the 16 
points). In this range, the current data is in complete disagreement with the predictions 
of the USJ model, with a significance of 99.8% according to the 2D K-S test. We checked 
the robustness of this result by varying our assumptions and by examining a few possible 
selection effects in the viewing angle 9, and found it to be robust. Thus, although the 
relatively small size and the inhomogeneity of the current sample prevent us from drawing 
a definite conclusion at this stage, we find that the current data disfavor the USJ model. 
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Fig. 1. — The 2D distribution density, n(z, ln9), of the GRB rate as a function of redshift 
z and viewing angle 9, as predicted by the universal structured jet (USJ) model. The white 
contour lines confine the minimal area which contains 1 o of the total probability. The circles 
denote the 16 bursts with known z and 9 from the sample of Bloom et al. (2003). (a): The 
parameters of the models are similar to these of PSF03. This figure is the 2D realization of 
their Fig. 1. (b) Here we use a limited range in redshift, 0.8 < z < 1.7 (containing 10 out of 
the 16 data points), in order to minimize redshift selections effects and reduce the sensitivity 
of the results to the unknown GRB rate. We take into account 20% measurement errors in 
\\\9 {(Ti n B — 0.2) and a lognormal distribution in T that we deduced from observations (Eq. 
5). 



